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PLOTTER  SUBROUTINE  FOR  METEOROLOGICAL 

FLIGHT  PATH  (PATH) 

Heinz  H.    Grote 

This  paper  describes  a  subroutine  PATH  that 
works  in  connection  with  the    plotter  subroutine 
CRTPLT  to  provide  a  two-dimensional  plot  of 
a  meteorological  flight  path  with  wind  vectors, 
time  notations,    and  a  legend. 

1.     INTRODUCTION 
The  subroutine  PATH  works  in  connection  with  the  plotter  subrou- 
tine CRTPLT  (Tucker,    19&8)  to  provide  a  two-d;mens  ional  plot  of  a 
meteorological  flight  path  with  wind  vectors,    time  notations,    and  a  legend 
Figure  1   shows  an  example  of  a  plot.      The  legend  provides  the  following 
reference  information: 

(1)  The  flight  number  LID  is  expressed  as  a  six-digit  number 
and  an  index  letter  according  to  the  convention     used  by 
ESSA' s  Research  Flight  Facility.      The  first  two  digits 
represent  the  year,    the  next  two  the  month,    and  the 

last  two  the  day  of  the  flight.      The  letter  stands  for  the 
particular  airplane. 

(2)  The  altitude  ALT  in   kilometers  is  a  reference  value 
brought  in  from  the  main  program. 


CENTER 

LONGITUDE  106.372  DEGREES 

i — f 

=  5  KM 

LATITUDE  -58.543  DEGREE? 

WIND1 ' 

=  10  N/S 

FLIGHT  N<3  69G815  B 
ALTITUDE   1.2  KM 

Figure  1.    Example  of  plot. 

(3)     Longitude  LOC  and  latitude  LAC  in  degrees  refer  to  the 
small  circle  in  the  middle  of  the  plot.      These  reference 
values    are  brought  in  from  the  main  program.      Together 
with  the  scale  notation,    every  point  on  the   plot  is  deter- 
mined by  this  location  reference  under  the  assumption 

that  the  earth's  surface  can  be  represented  by  a  plane. 

j 
The  relation  between  distances  in  geographic  degrees 

of  longitude  d\  and  latitude  dcp  to  linear  distance  values 
a  and  b  in  kilometers  is  given  by  the  equation 


a  =  C.    dX 

b  =  C  .    dcp  cos  cp, 
where  a  and  b  are  in  (km);   d  X,    dcp,    and  cp  are  in  geo- 
graphical degreees;  and  C  is  approximately  a  constant 
with  the  value 

C  =   111.  32 


geographical  degree 

(4)  The  reference  length  is  equal  to  a  distance  that  is  de- 
termined automatically  by  the  program  in  such  a  way 
that  the  following  scale  lengths  are  obtained  for  a 
corresponding  span  (SPAN     introduced  from  the  main 
program)  of  the  flight  path  in  its   largest  dimension 

Dimension  Scale  Length 

55  km  2    km 

137.5     km  5    km 

27  5  km  10    km 

above         27  5         km  20    km 

(5)  The  wind  reference  length  is  equal  to  a  wind  speed 
that  is  determined  automatically  by  the  program. 
Scale  lengths  for  a   maximum  wind  speed  (WM   intro- 
duced from  the  main  program)   are: 


Wind  Speed  Scale  Length 

10     m/sec  5    m/sec 

20     m/sec  10    m/sec 

50     m/sec  25    m/sec 

above     50     m/sec  50     m/sec 

(6)      A  wind  reference  arrow  is  provided  with  its  origin  at 

the  reference  coordinate  point.      Scaled  according  to  the 
wind  reference  length,    its  value  is  brought  in  from  the 
main  program  as  WXA,    WYA. 

If  this  vector  is  omitted,    the  represented  wind  vectors 
stand  for  the  wind  as  measured  at  each  flight  point.     If 
this  vector  has  a  value,    it  generally  represents  the  mean 
wind  over  the  plotted  area  or  a  closed  loop  withinthis  plot; 
the  wind  vectors  then  give  the  deviation  of  the  wind  at  each 
location  from  this  mean  wind. 

2.      VARIABLES 
The  variables,    listed  in   alphabetical  order,    that  enter  the  subrou- 
tine  are  as  follows: 

ALT  a  representative  flight  altitude  (km)  printed  under  legend  (2) 

with  one  digit  after  the  decimal  point,    so  that  hundreds  of 
meters    can  be  represented. 
IFS  a  six-digit  number  that  determines  the  scaling.      The  center 


two  digits  are  without  significance  and  serve  only  as  a  separa- 
tor.      The  first  two  digits  are  the  numbers  inserted  under 
legend  (4)  as  the  plotting  scale  values.      The  last  two    are 
the  numbers  inserted  under  legend  (5)    as  the  wind  scale 
values.      If  IFS  is  zero,    automatic  scaling  is  used,    as  ex- 
plained under  (5)  and  (6). 

IN  the  number  of  coordinate  points  X,    Y. 

INCW         determining  the  spacing  of  the  wind  vectors  on  the  plot. 

INCW  =    1  means  that  a  wind  vector  is  drawn  for  every 
coordinate  point. 

INCW  =  2    indicates  that  a  wind  vector  is  drawn  for  every 
second  coordinate  point. 

JT  the  time  array,    consisting  of  six  digits.      The  first  two 

stand  for  the  hours  (in  24-hour  notation),    the  second 
two    digits  for  the  minutes,    and  the  last  two  for  the 
seconds.      The  capacity  of  the  array  is  500. 

LID  the  identification  array,    containing  eight  places. 

LID(l)     consists  of  two  digits  and  represents  the  most 
significant  digits  of  the  year  of  the  flight. 

LID(2)    is  the  number  of  the  month  incremented  by  1000 
in  order  to  maintain  leading  zeroes  in  printing. 

LID(3)     is  the  number  of  the  day  incremented  by   1000 
in  order  to  maintain  leading  zeroes  in  printing. 

LID(4)    is  the  airplane  identifier  in  Hollerith  notation.      (The 
airplane  identifier  B  thus  has  to  be  expressed  as   1 
HB.  ) 


LID(5)  is  the  day  in  two-digit  presentation. 

LID(6)  is  the  first  eight  letters  of  the  month  in  Hollerith 
notation.      (The  month  of  August  is  thus  expressed 
as  6HAUGUST.  ) 

LID(7)  is  the  extension  of  the  month  from  the  ninth   letter 
in  two-digit  Hollerith  notation.     If  the  month  has 
less  than  nine  letters  LID(7)  =  2  H/\A.     (For  the 
month  of  September  LID(7)  =  2  HR/\  .  ) 

LID(8)  is  equal  to  LID(l)  and  contains  the  last  two  digits 
of  the  year. 

LAC  the  reference  longitude  in  geographical  degrees  as  explained 

under  legend  (3). 

LOC  the  reference  latitude  in  geographical  degrees  as  explained 

under  legend  (3). 

MT  the  spacing  of  the  records  in  seconds.      This  value  deter- 

mines the  time  lettering  on  the  plot.      The  first  time  on  the 
plot  is  a  four-digit  number  representing  hours  and  minutes. 
The  following  time  values  are  generally  expressed  by  only 
two  digits  representing  minutes.      Four  digits  are  printed 
only   at  intervals  of  MV/100.      Time  values  are  printed  only 
for  every  MV/100  minutes.      The  values  of  MU  and  MV  as 
they  correspond  to  MT  values  are  as  follows: 

120 
10000* 
2000 

^Represents   1  hour  due  to  the  conversion  from  the  60-bit 
time  system  to  the  decimal  system. 


MT 

5 

10 

20 

30 

60 

MU 

500 

1000 

2000 

3000 

10000* 

MV 

100 

100 

200 

500 

1000 

SPAN 


WM 


wx 


WXA 


WY 


WYA 


X 


the  largest  extension  of  the  flight  path  in  kilometers 
in  either  longitudinal  or  latitudinal  direction.      This 
value  is  used  for  the  automatic  scaling  as  explained 
under  legend  (5). 

the  largest  wind  speed  (m/sec).      This  value  is  used 
for  the  automatic  scaling  as  explained  under  legend  (6). 
an  array  of  500  words  that  contains  the  x- component  of 
the  presented  wind,    either  absolute  or  relative  to  the 
mean  wind  as  explained  under  legend  (7). 
the  x-component  of  the  reference  wind  as  explained 
under  legend  (7). 

an  array  of  500  words  that  contains  the  y- component 
of  the  presented  wind. 

the  y- component  of  the  reference  wind  as  explained 
under  legend  (7). 

is  an  array  of  500  words  that  contains  the  x-coordinates 
of  the  points  of  the  flight  path  in  a  plane  coordinate  sys- 
tem with  distances  from  the  reference  point  in  kilometers, 
an  array  of  500  words  that  contains  the  y- coordinates  of 
the  points  of  the  flight  path  in  a  plane  coordinate  system 
with  distances  from  the  reference  point  in  kilometers. 


3o      PROGRAM  TESTPATH 

To  test  the  subroutine  PATH,    the  program  TESTPATH  was  run 
through  the  ESSA-Boulder  computational  facility.      The  resulting  plot 
is   shown  in  figure   1,    and  the  essential  parts  of  the  printout  are  pro- 
duced in  appendix  C. 

The  variables  in  TESTPATH  and  PATH   are  compared  in  table  1 

Table  1.     Comparison  of  Variables 


TESTPATH    PATH 


TESTPATH 


PATH 


ALA 

LAC 

ALO 

LOC 

ALT 

ALT 

IFS 

IFS 

IW 

INCW 

JJ 

JT 

LT 

MT 

NM 

IN 

OX 

X 

owx 

wx 

OWY 

WY 

SPAN 

SPAN 

WXA 

WYA 

WMAX 

WM 

In  TESTPATH  ZDR  is  the  conversion  factor  to  express  degrees  in 

radians.      The  test  path  consists  of  three  circles  (loops   1  and  4,    2  and  5, 

3  and  6).      Loop  7  establishes  the  winds.      (Since  IW  =  2,    only  every  sec- 

cond  value  is  being  used  in  this  presentation.  )      WXA  and  WYA  are  the 
mean  values  established  by  dividing  the  sums  SWX  and  SWY  by  the 
number  of  points  (100).      The  second  OWX,    OWY  is  then  the  relative 
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wind  after  subtraction  of  the  mean  wind  WXA,    WYA. 

4.      FURTHER  DETAILS  ON  SUBROUTINE  PATH 
ID  the  identification  printed  on  the  microfilm.      It  contains 

the  name  of  the  user,    his  telephone  number,    and  the  pro- 
ject number  in  Hollerith  notation. 
UK  =  2J    .      Since  Hollerith  quantities    with  three  characters  are 

used  in  the  plot  routine,    integer  values  have  to  be  shifted  by- 
five  characters  for  conversion       (Control  Data  Corporation, 
1967,     1.1,    "Word  Structure"): 


666166666 


IMX  =        IN  +    1.      The  number  of  the  reference  point 
X(IMX),    Y(IMX)  -  0.  0. 

XD(1)  =     16    *  KK.    The  plot  contains   1024  points.      After  subtraction  of 
the  legend,    880  points  remain  in  the  y-direction.      In  the 
x-direction  some  space  is  needed  for  time  notations,    so  that 
a  square  of  880  points  is  used  as  a  basis.     Since  the  length 
of  the   reference  line  is  32  points,    the  span  that  can  be  re- 
presented is 


880 
1024 


*  32  *  KK 


KK  =  2 

KK  =  5 
KK  =    10 


SPAN  =  55  km 
SPAN  =  137.  5  km 
SPAN  =    275  km 


Because  of  the  legend  the  plot  is  not  distributed  symmetri- 
cally.     The  legend  takes  about  14  percent  of  one  dimension, 
so  that  the  upper  part  takes  (100-14):  2  =  43  percent  above 
zero^and  the  total  lower  part  extends  to  (100-14):  2  +   14  =  57 
percent  below  zero.     Since  the  total  span  is  32  *  KK,    the  upper 
y-limit  is  now  0.  43  *  32  =  13.  76  *  KK,    while  the  lower  limit 
is  now  0.  57  *  32  =  18.  24  *  KK.      Thereby: 

XD  (1)  =  16.  00  *  KK 
XD(2)=    -16.  00  *  KK 
YD  (1)  =  13.76  *  KK 
YD  (2)  =  -18.24  *  KK 

The  variables  used  for  plotting  are: 

X,  Y,  IN       the  coordinates  of  the  points  representing  the  flight  and 
their  number. 

3H+0L; 

X,  Y,  circle  marks     0    denoting  points  on  the  flight  path  and  in 

IMX  the  center. 


XD,  YD, 

ID  limits  of  the  plot  and  identification  of  the  customer. 

XL,  YL, 

KS,  KA        lettering  of  the  legend. 

3H+0L; 

XS,  YS         scale  end  marks  and  scale  lines. 

XT,  YT, 

JZ,  K  time  numerals. 

XW,  YW      wind  vector. 


10 


Further  variables  in  the  subroutine  PATH  are: 

SSC  a  lettering  scale  factor.      The  size  of  the  letters  in  plot-size 

is   12  points,    while  the  total  number  of  possible  points  is 
1024.      The  total  span  is  32  *  KK.      Thus  the  scale  value 
for  one  letter  is 

SSC  =       — ^ —      *  32  *  KK 

1024 

SSC  =    0.  375  *  KK 
LS  a  switch  that  makes  it  possible  to  print  four  digits  of  time 

for  the  first  labeled  time  value  JC. 
JC  the  first  time  value  to  be  printed.     It  has  to  be  an  integer 

of  MU. 
NS  a  switch  that  bypasses  the  selection  of  two  digits  for  the 

time  print.      It  is  activated  for  the  first  time  value  to  be 

printed. 
IZ  the  time  values  in  Hollerith  notation.      If  only  two  digits  are 

to  be  printed  the  last  two  digits  are  blanks:   3H  +   0  a   . 
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APPENDIX  B 
Time  Scale  for  PATH  Lettering 


MU  MV 
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500   100   1.1.1.1.1 
41  min  of  data 


[Sec] 
01  02  MT 

1.1.1.1.1.1.1.1.1.1.1   [1605]    5 
2  digits  every  minute   4  digits  every   5  min 


1600 
1000    100   1 


01  02  03  04         05 

1 1 1 1 1    [1610]   10 


1  hour,  23  min  of  data 


2  digits  every  minute  I  4  digits  every  10  min 


1600  02  04         06  08  10 

2000   200   1..1..1..1..1..1..1..1..1..1..1    [1620]    20 
2  hours, 46  min  of  data  |  2  digits  every  2  min   I  4  digits  every  20  min 

1600  05  10  15 

3000   500   1.1.1.1.1.1.1.1.1.1.1.1.1.1.1.1    [1630]    30 


4  hours , 10  min  of  data 


2  digits  every  5  min. 


4  digits  every  30  min 


1600                  10                 20  30 

10000   1000   1  ....  1  ....  1  ....  1  ....  1 1    [1700]    60 


8  hours,  20  min  of  data 


2  digits  every  10  min 


4  digits  every  hour 


1600  20  40  1700 

10000   2000   1  ....  1  ....  1  ....  1  ....  1    ...  1  ....  1  120 

16  hours,  40  min  of  data     2  digits  every  20  min.   4  digits  every  hour 


MT 

5 

10 

20 

30 

60 

120 

MU 

500 

1000 

2000 

3000 

10000 

10000 

MV  v 

100 

100 

200 

500 

1000 

2000 
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APPENDIX  C 

TESTPATH  Printout 

PROGRAM  TESTPATH 

0  TENSION  LTO(  P  )  »  JJ(  ^0  )  ,0X  (  ^00  )   ,0Y(  ^n,  ,0yx  (*00)  »0WY(  «S00  ) 
1  FOPMAT  (^X»l0r|7 .3  ) 

9  FORMAT  (5X*WXA  =  *F7.^,^X*WYA  =  *c7 .  ?  ,  ^X*  SWX  =  *FiO.^* 
15X*5WY  =  *Fl0,3) 

DATA  (ALA=-48.5^8)»(  ALO=106*372  )  »  (  ALT  =  !  •  *>  *4  )  ♦  (  I  FS  =  0  )  »  (  IW=7)  » 
1  ( L 10=69 » lCC8> 1015 >1HR» 15  »6HAUGUST*2HH   »6CM,(LT  =  60)» 
2 (NM=10  0) » (SPAN=6C.C) »(WMAX=1 5.0) » (ZDR=n #01 745329253) 

DO  1   1=1 i 2 

JJ(1 )=(2257+I )*100 
1  CONTINUE 

DO  2   1=3,6? 

JJ( I )=(2297+I )*100 
?    CONTINUF 

no  3   I=63»l00 

JJ(  I  )s(_63  +  -I  )*100 

3  CONTINUE 

DO  4   1=1,27 

A  =  { 758-6*1 )*ZDR 

OX( I )=18.0*C0S( AJ-6.0 

OY(  I ) =1 P.^*SIN( A  1  +  17.0 
u    CONTINUE 

DO  5   1=28,73 

A=( ?07-4*I )*ZDR 

0X( I )=30.0*COS( A )-6.^ 

OY(  !  )  =-aO.O*SIN(  A  ) 

5  CONTINUE 

DO  6   I=74»l00 

A= (708-6*1 )*ZDR 

OX( I ) =18.0*COS( A  5-6.0 

OY( I )=18.0*?IN( A) -12.0 

6  CONTINUE 
SWX=SWY=0.0 

DO  7   j=l  ,50 

0WX( J)=1 0.0*COS( ( ?80-4*J)*ZOR )-l ?.o 
0WY( J )  =  15.0*5  IN ( ( 4  70-8*J)*ZDR )+  8.0 
SWX=SWX+OWX (J ) 

swy=swy+owy( j) 
7  continue 

pPTMT  1 ,OWX ,OWY 
WXA=SWX/50.0 
WYA=SWY/*0.0 
no  p   j  = l  ,50 
0WX( J)=OWX( J)-WXA 
0WY( J) =0WY( JJ-WYA 
8  CONTINUE 

PRI NT  ? , WX A , WY A , SWX , SWY 
P°INT  l iOWX ,OWY 

CALL  PATH  (LID»ALT , IW»NM» JJ»LT»SPAN ♦ 
1  W^AX,WXA,WYA,OX,OY  »OWX »0WY»AL0 • ALA » I^S ) 
ENH 
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APPENDIX  D 
PATH  Printout 

SUBROUTINE  PATH(LID»ALT» I  NOV* IN»JT»MT» 

1 S»AN  » WM  » WX A  » WYA  »  X  ,  Y  ,  WX  t'WY  *  LOC  9  LAC  9  I F S ) 
TYPF  REAL  LOOLAC 
DIMFNSION  IZ (4) »JT(R00) >KS(116)  »LID ( p  )  »MUT (ft )  »vvT(6)  >WX< 5  00)  9 

1  W  Y  (  q  0 0  )  ♦  X  (  r  C  C  )  ,  X  n  (  ->  )   9  XL  (  1  1  A  )  *  XS  (  ?  )   ,  X T  (  l  )  >  XV  (  7  )   *  Y  (  e  0 0  )  9  Y n  (  ?  5  , 

2  YL  (  1 16  )'  »  YS  (  ?  )  #  YT  (  4  )  t  YW  (  2  ) 

DATA  (  ID=33HH.GROTE    EXT. 6270    PN  21F503190) 

DATA  ( IJK=1073741824) 

DATA  KS=3H+0C»^H+0E 9 3H+0N j^H+OT ,3H+OF , 3H  +  OR 9 


l'3H-t-0F9^1-l  +  nL 
23H+0  93H+0 
^m+OL  9-SH  +  OO 

4^m+o  »o,o»o 

r^h+OD »3H+0F 
ft-3u+0=  ,'H+O 
7-^M-fCA  »^H+  Ol 
8^H  +  0  9^H+0 
P3H  +  DL.  »^>H  +  0A 
1 ? (3H+ 


»3H+0I  »3H+0G»3H+0H »3H+0T  9  3H+0  >3H+0N »3H+0O » 
»0,0,0,C,0»3H+0  t0» 

^H  +  OM^H  +  ^n  ,^m+0I  »^M+OT  ,^u+OU  ,^u+pn,3H+0F  9 

93H+0. ,0,^,0,3H+n  , 

,^H4-0r,^H  +  ^P,?H+0'r93M+0F,'3H  +  0c, 

jO^C^H+O  »3M+Q*»3H+0M» 

,3H4-CT  9  3H+0T  *3H+0T*3H+0U  *3H+0n  j^M  +  n^, 

9O  jOf^H+0,  9O  »3H+0  ,3H+0K  *3H+0M  , 

9  3H+0T9  3H+0I  9  3w+0T9  3H  +  OU9  3H-t-Cn,?H+OE9 


)  ,n,o, ^M+' 


»o»n»c-»?H+o  » 


^^M+Qn  ,?H+0r»3H+0G»  3H+0R  »^M+0fr  »3m+0F  *^H-».0S  » 
33H+0W»3H+0I  »3H+ON»3H+OD»6(3H+0     )  »3H+0=  9 
a^m+0  »0»0»?H+0  >3H+0M»^H+0/»3M+0S) 
HaTa  (MUTsT  00  ti  00  »?0 Ot «? 00  »i  000  » 9  000)  t.(MVT  =  «;00»l.  000*20-00  t 

1  •20^'n,1Or'r»n9lr'00ri) 

IMX=IN+1 
X(TV!X)=0.0 
Y  (  t  m  X  )  =  0  .  0 

IF  (  IFS.NP.O)   o7*9? 
07  Kf=IpS/l0000 
K1=KK/10 
K2=KK-K1*10 
L=tfs-'<:k.*i.OOOO 
L1=L/10 
L2=L-L1*1C 
GOTO  99 
Kl=2 
K2  =  0 

IF  (SPAN  .LE.  275.0) 
Kl  =  l 

Ic  (SPAN!  .L^.  137.  5) 
<1  =0 
K2  =  5 

IF  (SPAM  .LE. 
Ll=5 
L2=0 

IF  (WM.L'F.50)  49  7 
LI  =? 
L?  =  5 
IF  (VM.LF.70)  5  97 


98 

1 
? 


1  93 

2*3 
5  5.0)   K  2  =  2 
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6*7 


5  Ll=l 

L2  =  0 

IF  (WM.LE.10) 
f  L1=0 

L?  =  ^ 
"7  K<s1  0*^1  +  v  ? 

99  XD(1)=  16.00*KK 

xo(?)=-i6.no*<< 

YD(1)=  13.76*KK 

Yn'(  ?  )  =-1  R  ,?4#KK 

SSC=0.375*KK 

CALL  CRTPLT  ( Xp  »  YP , A » ID , 2 ) 

CALL  rPToi_T  (0»0»*H+0{»4»5) 

CALL  CPTDLT  (X«YfIMX»l*l) 

XS  M  )  slO.5*^^ 

X$(2)=11.5*KK 

YS(1 )=YS(2)=-14.8P*KK 

TALL  TPTPLT  (  C,C»-3H  +  0L»4  »  «  ) 

CALL  CRTPLT  ( XS » YS »2 » 1 » 1 ) 

CALL  CRTPLT  ( XS  »  YS * ?  » n » 9 ) 

YS (1 )=YS(2)=-15.63*KK 


(XS»YS»2tl»l) 
(XS*YS»2  ,0,9) 


CALL  TPTPLT 
CALL  CRTPLT 
CALL  CRTPLT 
JM=IN/INCW 

JF=JM+l 

JP=JF*INCW 

X(  JF)=Y(  JF)=CO 

WX( JC)=WXfi 

WY( JF)=WYA 

DO  8   J=1,JF 

JI=J*INCW 

XWfl  )=X( J  I  ) 

XW(2)=X( JI )+WX( J)*KK/L 

YWd  }=Y{  JI) 

YW(2) =Y< J I )+WY( J)*KK/L 

C&LL  CRTPLT  ( XW»YW %?.  »0»9  1 

CONTINUE 

CALL  CRTPLT  (0*n»OfO»8) 

r.^LL  rRTPLT  ( X »Y» I N» 0»1 J 

SSC=0.375*KK 

ASSIGN  82  TO  LS 

ASSIGN  85  2  TO  NS 

JTsMT/lO+i 

T  p ( V  T    .LT.       iO)        tT=i 

IF(MT    .GP".       60)        IT  =  S 

IF(MT     .GE.     120)        IT=6 

MU=^UT( IT ) 

MV=MVT (IT) 

JC=CJT(1  )-l  )/MU*MU  +  MU 

po    RQ        1=1 t IN 
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GOTO    LS»(82»84) 
82     IF     (  JT(  I)     .GE.    JO        83»89 
R-a     ASSIGN    84    TO    LS 

GOTO    85 
84    IP(MOD(  JTU  )  »MU)     .PQ.    0)       R^>89 
8*     JX=JT( T)/]0000 

JY=( JT( I )-JX*10000)/lC0 
I<A=JX/10 
IKB  =  JX-KA*10 
I<C=JY/in 
IKD=JY-IKC*10 
IZ(1  }=IKA*IJK 
lZm  =  IKB*IJK 
IZ(3)=!KC*IJK 
IZ(4)  =  IKD*UK 
GOTO    NSt(851»852") 
SM     jP(MOn( JT( ! ) »MV)     .FQ.    0)        87,86 
852    ASSIGN    851    TO    NS 
GOTO    87 
s*     T  Z ( 1  )  =  T  Z ( *) 
IZ(2)=IZ(4) 
IZ(3 )=*H+0 
TZ(4)=^+C 

87  XT(1)-X( I )+SSC 
XT(2)=XT(1 )+SSr 
XT(3)=XT( ?)+SSC 
XT(4)=XT(3)+SSC 

YT(1  )=YT(7)=YT<3)=YT(4)=Y( I )-0.!5*SSC 

DO    88       K=l»4 

TALL    rRTPLT     (0,O,lZ(<) »1»51 

tail    roToLT     (XT(K)»YT    U)»1»0»1) 

88  CONTINUE 
po    CON'TTNUP 

Do    15         KA=1*6 
]  k    XLdi<A)  =  (<&-4)*SSr 
DO    16  |<A  =  7»24 

16  XL(<A)=(KA-48)*SSC 
DO    17  KA=25»49 

17  XL  (K  &  )  =  (  KA-^>7)^SSr 
DO    18         KA=50,56 

T8    XL_(KM  =  (K&-16)*SSC 
DO    19  KA=57,73 

19  XL(KA)=(KA-98)*SSC 
DO    20         KA=74»98 

20  XL(KA)=(KA-86)*SSC 
DO    21  KA=99,116 

2\     XL(KA)=(<A-75)*SSC 
DO    22         KA=1»6 

22  YL(KA)=-38*S^C 
DO    2  3         KA=7»56 

23  YL(KA)=-40*SSC 

r>Q    24  < A  =  57,  116 
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?Ia     Yl  (  K&  )  =-£.?*5SG 

KS(17)=LID( 1)/10*IJK 

IF     (KS(17).EQ.O)     KS(17)=1H 

KS<18)=(LlD(l )-LlD(l)/lO*lO)*IJK 

NM=LID(2)-LID( 2 )/ 100*100 

icS  (  i  °  1  sNM/1  C*I  J1*' 

KS( ?0)=(NM-NV/]0*10)*IJK 

NN  =  LID(3)-I_ID(3)/100*1C0 

<S(  ?•)  )  =  NJ  N  /  T  0  *  T  J  W 

KS(?2)  =  (NN-NN/lO*l.  0)*IJK 

KS(24)=LlD(4)/4096 

IF     (LOG     .LT.     0.0)        <c(?4)=?H+0- 

l_OC  =  APS(LOC  ) 

[_0A  =  TNT  (  L^r  ) 

L01=LOA/lCO 

KS(35)=L01*IJK 

Ic  (L01.FQ.O)   KS(35)=?H+0 

LOR=LOA-L01*lOO 

LO?=LOR/10 

KS(36)=L02*UK 

IF  ( (L01  .EQ.  0) 

LOD=L03-L02*10 

KM37)=L00*IJK 

LOF=IMT(  1 000.0*LOr  ) 

LOF=LOF-(LOF/1000)*1QOO 

LOG=LOF/100 

KS(39)=LOG*IJK 

LOH=LOF-LOG*100 

LOI=LOH./lO 

KS(40)=LOI*IJK 

LOj=L0u-Lni*1 0 

KS(41)=L0J*IJK 

<S(52 )=K1*IJK 

IF  (Kl.FO.O)  KS(5?)=3H+0 

KS(53)-K?*UK 

jalt=int( 1c*alt) 
ju=jalt/ioo 

<S(6"M=JU*TJK 

IF  (JU.EO.O)  KS(67)=3H+0 

JA=JALT-100*JU 

JV=JA/10 

KM68)=JV*IJK 

JW=JA-JV*10 

KS(70)=JW*IJK 

IF  (LAr  ,LT.  0.0) 

LAr=A°^(L*C) 

LAA=INT(LAC) 

LA1=LAA/10 

KS(85)=LA1*IJK 

IF  (LA1.FO.0)   <S(85)=2H+0 

L^c'  =  i_/\A-Lfl1*l0 

KS(86>=LAB*IJK 


AND.  (L02  .EG.  0))   KS(36)=2H+G 


<5(  PL  )  =^M  +  0- 
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Lac=INT(  1.000. 0*LAC) 

LAr  =  LAc--(LAr/T00  0)*i000 

LAG=lac/100 

<S(88)  =LAG#I  JK 

LAH=LAF-LAG*100 

LAI=LAH/10 

KS(89)=LAI*IJK 

LAJ=LAH-LAI*10 

KS(90)=LAJ*IJK 

KS( 11! )=L1*IJK 

If7  (L1.EQ.0)  KS(lll)=3H  +  0 

<S( 11?)=L2*IJK 

On  25    KA  =  1  ,116 


_  7      x  J.  V.J 

CALL  CRTPLT  (  0  »0»KS'(  KA)  »  1  »5  ) 

(XL(KA)  »YL(KA) ,1»1»1 ) 


CALL  CRTPLT 
CONTINUE 
CALL  CRTPLT 
RETURN 
END 


(0,0,0*0,20) 


GPO  855-  538 
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